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Abstract 

The full evolution during the late periods of the condensation pro- 
cess is described in the analytical form. The process is split into several 
periods and for every period the simple approximate solution is given. 



1 Introduction 

The first order phase transition is characterized by the temporal duration - 
the process of the condensation lasts in time and consists of several charac- 
teristic periods. The kinetics of the first order phase transition is one of the 
actual problems in the phase transformations. 

Ordinary it is supposed that the final stage (period) of the phase tran- 
sition is the stage of coalescence. This name goes from the theory given by 
Lifshic and Slezov pp. This terminology is not absolutely correct - formally 
the coalescence means the adhesion of the embryos. But the Lifshic-Slezov 
(LS) consideration does not take into account the adhesion, the evolution in 
the LS picture is the competition between the embryos through the exhaus- 
tion of the vapor environment. One has to clarify this difference. In the case 
of adhesion we shall speak about the coagulation and when the Lifshic-Slezov 
mechanism takes place we shall speak about the over-condensation. 

The over- condensation means the competition between the already formed 
embryos when the embryos of the relatively big sizes eat the embryos of 
the relatively small sizes not directly but through the exhaustion of the 
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metastable old phase until such a degree that the mentioned small embryos 
become the pre-critical ones and begin to eject the molecules in the vapor. 
These embryos begin, thus, to disappear. 

Beside this opportunity it is possible to see in the open systems the situ- 
ation when all embryos continue to grow until all volume of the system will 
be transformed into a new phase. Certainly this situation takes place when 
there is an effective source of the substance in the old phase. This is also an 
alternative possibility for the final of the phase transition. So, we see that 
the final stages of the phase transformation can be very different. 

Here we shall speak about the over-condensation and consider the closed 
system. The effect of coagulation is not taken here into account. The tem- 
perature of the system is supposed to be fixed. The mentioned restrictions 
are not the crucial ones. 

The methodology of analysis of such systems is given by Lifshic and Slezov 
in pp. Strictly speaking the results of Lifshic and Slezov can not be directly 
applied for the systems with the diffusion regime of the substance exchange 
(as it was done in [T]) because in this case around every embryo there will 
be a profile of the metastable (old) phase substance. The profiles around 
different embryos overlap. This changes the rates of growth. Then one has 
to consider an interaction between the density profiles. So, one can not 
directly take the integral form of the balance equation as it is done in the 
theory of Lifshic-Slezov. The last task to describe the interaction of profiles 
is extremely difficult and the precise analytical solution of this problem is 
hardly possible. 

The ignorance of the diffusion profiles leads to the necessity to consider 
the situation with the free molecular regime of interaction when there is no 
profiles. This was done by Wagner [2] and we shall follow this regime of the 
substance exchange below. 

The fundamental fact is that the LS asymptotics was confirmed exper- 
imentally (see [3], [I], [5]). Here one has to stress that according to the 
notation made in [6] the accuracy is not very high - one can speak only 
about the qualitative confirmation of the form of the spectrum and about 
the precise confirmation of the power-like law of evolution of parameters of 
the spectrum (in some sense the power-like law is rather evident and it can 
be derived from some simple qualitative remarks). Also it is necessary to 
stress that according to the references from j6] the form of the size spectrum 
is more sharp than it follows from the leading term of the LS theory. 

The next step in modification of the LS approach was done by Osipov and 
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Kukushkin in [B] where the "regular asymptotic" was constructed. Since the 
first two terms in asymptotic decomposition in the LS theory are universal 
(independent from the initial size spectrum) it is possible to introduce such 
variables (coordinates) where already the initial approximation is equivalent 
to the first two terms of the LS decomposition. So, the modification of Osipov 
and Kukushkin seems to be natural. But this modification is important 
because it allows to speak on the level of the regular asymptotics already 
from the very beginning. 

The LS asymptotics has an amazing feature - it allows to establish the uni- 
versal form of the size spectrum. One can admit that the power of metasta- 
bility has some universal asymptotics, but one can hardly believe that the 
spectrum of the embryos sizes is universal. This will be the matter of dis- 
cussion is this paper. Below it will be shown how the real spectrum of sizes 
approaches the form given by LS theory. 

2 Some remarks initiated by technique of Lif- 
shic and Slezov 

Now we shall consider the simplest model leading to the over-condensation. 
At first we shall consider the embryos with linear sizes which strongly exceed 
the critical size. We shall call them as the supercritical embryos. The growth 
of the supercritical embryo in time t is given by the following equation 

dp = c 

dt ti 

Here p is the cubic root of the number of molecules inside the embryo (it plays 
the role of the linear size), ( is the supersaturation of the metastable phase, 
t\ is some characteristic time (this is simply the constant coefficient). The 
supersaturation is defined as the ratio of the surplus density of the metastable 
phase (with respect to the saturated phase) to the density of the saturated 
phase or the ratio of the real density of metastable phase to the density of 
the saturated phase minus one. 

The simplest asymptotic correction to the asymptotic law of growth will 
be the following 

^=i(i_„-i) 

at t\ 
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Here 

Pc 

is the ratio of p to the critical size p c to have the zero rate of growth for 
the critical embryo. It is quite possible that this ratio has been changed 
by an analogous ratio in some power but this does not lead to the essential 
reconsideration of this approach. 

One can give the interpretation of the last correction term as the lead- 
ing term in the asymptotic expansion in inverse powers of the linear size. 
It is possible that this decomposition starts from the non-integer power as 
it takes place in construction of series for solutions of the linear second or- 
der differential equations. It is also possible that the law of growth for the 
supercritical embryos has some power of the linear size which grows with a 
constant velocity. Namely this situation takes place in the diffusion regime of 
the metastable phase exchange. In both cases the technique described below 
can be applied. 

So, since one can see that 

<*>/*> _ = (i _ ,/„) 



(dp/dt) 



asymp 



or there is some function of u in the rhs it is convenient to consider du/dt 
instead of dp/dt. For the derivative du/dt we get an additional term 



dp du dp c 
dt ^ c dt dt 



which is linear on u. So, 



du C/ n -u d Pc 

— = —{l—U ) — U — ; — 

dt ti 1 dt 

In the LS theory there exists a hidden supposition that p c depends on time 
droningly. This supposition was reexamined in approach which predicts the 
periodic formation of the tail of the size spectrum and then the consumption 
of this tail. But nevertheless this supposition is rather natural at least at the 
asymptotics. If we adopt this supposition then all values depending on time 
are the values depending on the critical size. At least one can consider the 
intervals of the monotonous dependence. 
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The generalization of this law on other regimes of growth leads to 
on 

= ^(t)u a {l-u- l )-uUt) 

at 

The functions <fii, <p2 are the functions of time. The function <pi can be 
excluded by transition to the rescaled time r. The same can be done with 
u a . 

The right hand side of the last equation (let it be T) can not attain 
asymptotically the positive value at some argument - then the substance 
balance will be violated. The balance equation will be violated also when 
the rhs is negative at all values of argument u. Hence, the unique possibility 
is to touch the zero level at the main maximum. Then it is necessary that the 
rhs as a function of u has a maximum and then by the correct choice of ip one 
can put the maximum to the zero level. The class of functional dependencies 
which allow this operation forms the class of dependencies where the LS 
technique can be applied at least formally. The fine unjustified supposition 
here is the monotonous approach (in time) of the maximum to the zero level. 
It is quite possible that this maximum oscillates near zero - at some moments 
is it greater than zero, at some times it is less than zero. This possibility is 
considered in |8j. Here we do not consider the oscillating scenarios. 

For the function ip we get 

<P = 7o 

where 70 is some constant. Even when this supposition is adopted the velocity 
at the maximum is zero and the balance equation is violated. Then it is 
necessary to be: 

¥? = 7o(l + £ 2 (r)) 

with e at r — >■ 00. 

The reason that 70 is a constant lies in the form of equation for du/dr. 
It would be interesting to consider equations which do not belong to the 
established form and to get the evolution analogous to the over-condensation. 

Let the argument of the main maximum be u m . One can see that at 
ip = 70 the embryos with u > u m can be dissolved only at the infinite time 
(the diffusion is not considered here). 

Now we shall establish the function e{r). To find this function one has 
to analyze the solution near the maximum of V. We introduce the variable 

z = (u- u m )/e 
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Then we see that 



du/dr = d(u — u m )/dr = dze/dr 



edz/dr + zde/dr 



and 



dz 



r(7o(i + ^ 2 )) 



zk\ 




where k\ is the coefficient equal to de/dr. In the last equation it is necessary 
to express u in the function V through z and e. Then we rescale the time to 
exclude the coefficient in the rhs and then we get an equation which allows 
(or does not allow according to the reasons mentioned above) the analysis 
analogous to the already made one. If equation allows such analysis one can 
get the correction term of the asymptotic. 

One has to get the equation binding de/dr and e. In the power-like laws 
of growth this equation is trivial and has the form f(e) de/dr = const with 
some known function /. It is interesting whether one can get any more 
complex form of this equation. It is quite possible that some more complex 
form can lead to some new physical effects. Here we do not consider this 
question and assume that the dependence e on r for the correction term is 
established. 

Having calculated one correction term after another we establish all asymp- 
totic series until the step when we can not perform the procedure described 
above because the obtained equation does not allow to put maximum to the 
zero level. In the natural physical situations only the initial and the first 
correction term can be established. It is worth seeking the situations where 
one can make more steps. 

In the Osipov-Kukushkin approach we get the correction term already 
as the initial approximation since one can choose such variables where the 
correction term can be treated as the initial one. Certainly, it is possible 
to continue this procedure but it appears that under the natural regimes of 
the substance exchange it is impossible to choose parameters in equations 
for correction terms that the main maximum will touch the zero level. It is 
interesting to find situations where all corrections are the universal ones. 

One can see that the account of the diffusional term in the evolution 
equation can not lead to another asymptotics because the diffusion process 
along u becomes negligible at the infinitely big time. In correction terms it is 
necessary to check that the diffusion process is really negligible in account of 
corrections. It would be interesting to seek the situations where the account 
of the diffusion process changes the correction terms. 
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The unpleasant fact for the LS theory is the following one. In the zero 
approximation the time of the dissolution of the embryo with u > u rn is 
infinite and has to be refined. But in the first approximation the time of 
dissolution of every finite embryo is the finite one and it means that every 
finite size spectrum will be dissolved (dissolved at the finite time). After the 
spectrum is dissolved the balance equation will be certainly violated. 

The mentioned difficulty has a fundamental character which is confirmed 
by further constructions in frames of the LS theory. It is reasonable to 
connect the finite time of dissolution with the possibility to find the universal 
form of the size spectrum. The authors of the LS approach treat the universal 
distribution as the distribution which is the limit of the relaxation process. 
But it evident that namely the initial form of the size spectrum determines 
the whole further evolution. Under the regular and only under the regular 
law of growth we have 

p(p,t)dp = f(p',t')dp 

where p and / are the old spectrum and the new spectrum of sizes and p at 
the moment t has to come by the regular growth into p' at t'. Then 

i<M>g = /W,o£l«w 

Since for dp/dt we have the concrete given expression it is hardly possible to 
change the form of the size spectrum to come to the universal form of the 
LS theory. If one determines the velocity of growth from comparison of the 
initial form of spectrum and the final universal form of spectrum it leads to 
the unpleasant contradiction. 

The alternative is the following: to have the size spectrum prescribed by 
the initial distribution or to see the leading role of diffusion (at least at some 
time). 

It appears that the universal form of the size spectrum has absolutely 
another sense - this form of spectrum is such a form which corresponds to 
the already established asymptotics for the critical size precisely or at least 
ensures the optimal relaxation to the established asymptotics for the velocity 
of growth. With the real form of the size spectrum this asymptotic form of 
the size spectrum has no direct connection. Then the asymptotic for the 
velocity of growth corresponding to the ideal size spectrum will be never 
attained at finite time. 
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3 Asymptotic form of the size spectrum 



One can get the universal form of the size spectrum in the LS formalism 
rather simply. In the LS approach it is supposed that the main quantity of 
the substance is in the region u < u, m . Then the balance of the substance 
leads to 



1-UmPc 

/ p f(p,t)dp = const 
Jo 



where f(p, t) is the distribution function. From this function we come to the 
distribution over u, namely to <f>(u,r). Then 

(fi(u,r)du = f(p,t)dp 

The substance balance can be rewritten as 

p c / u <f)(u, r)du = const 
Jo 

The form of the function can be determined from the continuity equation 

d(f> d(pv u 
dr du 

where v u is the velocity of the growth for u. The solution of this equation is 
rather simple 

<p{u,T) = 6(t - t{u))/v u 

where t(u) is the time for the embryos to attain u. Instead of the last function 
one can write the function of u — u(t), where u{r) is the size attained at the 
given time. The sense of solution is that the dependence over two variables 
is performed through the dependence over one variable. The concrete form 
of this dependence is determined by initial conditions (in the LS approach it 
is detemined by asymptotic relations). 

The transformation to the variable u is made to ensure the constant value 
of the upper limit of integration. 

Ordinary the concrete form of the functional dependence of 9 has to be 
determined from the initial conditions. But in the LS formalism this form is 
determined from the asymptotic balance equation which can be written as 



p / u 9(t — r(u))/v u du = const 
Jo 
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The function v u depends only on u, but not on r. This is the consequence 
of the fact that in the zero approximation for du/dr one can take 70 which 
excludes the dependence on time. 

The function p\ at the already established asymptotic is the known func- 
tion of time. This dependence has to be cancelled by the true choice of the 
functional dependence for 9. Namely this cancellation is the recipe to choose 
the true form of the function 9. 

The simple dependence (but may be not the unique one) is the following 

one 

9(t - t(u)) = g(exp(r - r(«))) 
where the function g satisfies the relation 

g(ab) = g l (a)g 2 {b) 

for any a b. As this function one can take the power function. Since we have 
to cancel only the leading term in asymptotic the choice of the power-like 
function is quite satisfactory. 

So, the solution is announced. But is it well justified? 

Really, the ordinary solution of the problem (not the asymptotic one) is 
very simple it is already presented in the form of the function 9, i.e. in the 
fact that 9 is the function of only one variable. Certainly, this function has 
to coincide with the initial spectrum of sizes at the initial moment of time. 
It is sufficient to calculate the origin of the given embryo taking into account 
the known supersaturation, then to take the initial distribution and divide it 
on the velocity of growth in the given point (or on the ratio of the velocities 
of growth in corresponding points). 

We shall speak here about the evolution scenario. This scenario is in 
contradiction with the LS scenario. It conserves the explicit dependence 
on the initial distribution. It is clear that the evolution scenario is more 
correct in comparison with the LS scenario. The LS scenario has at least two 
disadvantages 

• -The absence of the spectrum at u > u m 

• -The spectrum is determined to ensure the asymptotic, but it is more 
reasonable to get the asymptotic on the base of spectrum. 

The first disadvantage can be ignored by notation that every finite spec- 
trum sooner or later will leave the region u > u m . 
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The second disadvantage is more serious because it destroys the method- 
ology of the LS approach. It seems that it is the crucial point and the real 
size spectrum will differ from the result of the LS theory. 

But one can see the amazing fact - the similarity between the form of the 
LS spectrum and the results of experiment. We shall explain this similarity 
below. 

4 Approximate similarity between the real 
size spectrum and the results of the LS ap- 
proach 

We start to consider the strange fact - why the form of the size spectrum 
prescribed by initial conditions will resemble the universal result of the LS 
theory? 

At first one has to stress the low accuracy in the experimental determina- 
tion of the size spectrum form. As an example one can consider the style of 
discussion in j6] . In j6] one can see not only the account of the initial approx- 
imation but already of the first correction in the form of the size spectrum. 
The form of the size spectrum with the first correction essentially differs from 
the form in the zero approximation. This fact is very important and leads 
at least to two essential conclusions 

• - the correction term in the asymptotic expansion is not small at least 
for the size spectrum - one can not see the parameter with a property: 
the small value of this parameter leads to the size spectrum in the zero 
approximation. Moreover, one can state that there is no such parameter 
because the size spectrum in the first approximation is universal one. 
This shows that there is no reliable way to determine the size spectrum 
because the next correction can change the form of the size spectrum 
radically. 

One can also give an interpretation which is not favorable for the LS ap- 
proach - the modified zero (Osipov-Kukushkin) approximation is only 
the starting point and all further terms of decomposition depend on 
the initial size spectrum. So, there is no reason to speak about the 
universality. 
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• - the accuracy of the experimental results is rather low. Really, the 
difference between the size spectrum in the zero (LS theory) and in the 
first approximation is very essential. Earlier the experimental results 
confirmed the LS spectrum in the zero approximation. Now experi- 
ments confirm the Osipov-Kukushkin result. It can be only when the 
accuracy is low. So, one can speak only about the experimental obser- 
vation of some tendencies in the form of the size spectrum. 

In the analysis of experimental results one has to note that some au- 
thors speak about the diffusional blurring of the size spectrum which 
corresponds to experiment. 

But one can show analytically that it is possible to neglect the diffu- 
sional term in the LS technique. It does not mean that in the evolution 
scenario one can neglect the diffusion - at least there are some periods 
when the diffusion is the driving force of evolution. Also one can see 
that in the case of the finite size spectrum one has to include diffusion. 

As the result of these facts one has to conclude that there exists a tail of 
the size spectrum (it is the exponential one - this is explained by the diffusion 
process) at u > u m (in any other region it can not be noticed). This tail is 
observed experimentally. What is the reason for this tail? The LS formalism 
can not give an answer. Below this answer will be given. 

Kinetics of the new phase embryos formation has some characteristic 
features which help to determine the characteristic features of the over- 
condensation. Roughly speaking, the "initial" size spectrum belongs to a 
narrow class of functional dependencies. So, we use the evolution approach 
and construct the sequence of stages for the process of condensation and 
over- condensat ion . 

The process of condensation (the periods before the over-condensation) is 
investigated in [7]. This investigation gave the total number of the embryos 
and the form of the size spectrum. For the further regular evolution it is 
necessary to know the first three (and the zero) momenta of the size spectrum 



or of the distribution function f(x) of the variable x, defined as the deviation 
of the coordinate p from the size z of the maximal value of this variable cor- 
responding to the embryo formed at the very beginning of the condensation 
process. 




11 



4.1 The period of the initial relaxation 

The relatively intensive formation of droplets stops at the relative decrease 
of the supersaturation equal to the value reciprocal to the quantity of the 
molecules in the critical embryo. 

The number of molecules in the critical value is a big value. The relatively 
small decrease of the supersaturation stops the nucleation (formation of new 
embryos) and later the size spectrum moves along the p-axis without any 
change of the form. To ensure the essential exhaustion of the metastable 
phase the spectrum has to move along p-axis for a rather long distance. So, 
at the end of this evolution the spectrum of sizes can be considered as the 
monodisperse one. 

The balance equation can be written in the following form 

$ = c + E ^JK-/ 3 -^i 

or 

which is the ordinary first order differential equation without the explicit 
dependence on the argument. So, it can be easily integrated. The result is 
the monodisperse spectrum and the relaxation of the spectrum coordinate 
to the critical size (or more correct the relaxation of the critical size to the 
coordinate of the spectrum). 

Here we use the law of growth for the supercritical embryos. It can be 
replaced by the precise law of embryos growth. Really, in the law 

dz = _( _ z3(, 
dt ti 2a 1 

we replace ( by $ — X^=o Ji^M ^^^i an< ^ S e ^ equation 

dt *i 2a ' 

which can be easily integrated. Here it is ignored that during the evolution 
the momenta p« will be changed which is considered as a correction. In any 
case we need only the initial approach to the critical size where the rough 
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monodisperse approximation is sufficient. Then there will be no problem 
with changing momenta. 

One can show that the account of the diffusional term will be essential 
only when the coordinate of the size spectrum is rather close to the critical 
size. This is the end of the relaxation stage and the beginning of the new 
period. 

4.2 The period of the diffusional blurring of the size 
spectrum 

The result of the previous period is the relaxation of the spectrum coordinate 
to the critical size. The spectrum resembles the delta-like function, The 
spectrum width 5p is many times less than the coordinate p or p c which is 
the spectrum coordinate. 

If there would be no diffusion then the spectrum will remain near the 
critical coordinate until the end of the whole evolution. But after the time 
of relaxation at the previous period the diffusion becomes the main driving 
force of the process. 

Kinetics of the diffusion blurring is rather simple and it is described in 
[S] • It is possible to approximate evolution by diffusion blurring without any 
regular growth with a boundary condition 

f(p = 0) = 

and the initial condition 

f(t = tiniuai) ~ S(p- p c ) 

The method to solve this problem is the combination of the Green func- 
tions at the infinite interval. The method of images allows to construct 
solution by addition of the negative gaussian in the symmetrical point. 

So, we write the diffusional equation in the following form 

^ = w*l 

dt dp 2 

Here W = W + + W~ is the generalized kinetic coefficient equal to the 
weighted sum of the adsorption coefficient W + and the ejection coefficient 
W~. One can take these coefficients in the critical point marked by the index 
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c. Then one can assume that W c = 2W C + . Taking into account the evident 
relation 

W + = Sv t /4 = 3p 2 (( + l)/t 1 

where v t is the mean thermal velocity, S is the surface square of the embryo 
one can determine the dependence of W c on p. Then it is necessary to go 
from p to a new variable r which is p in some constant power, i.e. p const . This 
transition excludes the dependence of the diffusion coefficient on the size at 
least asymptotically. It occurs at pdp ~ ds = dp 2 . 

The Green function at the infinite interval is written in the following form 

G( S ,t| So ,t )~exp(-^2l_) 

where D is the diffusion coefficient (this is the known constant value), s is 
the point of appearance of elementary disturbance at the moment t , s is the 
point of observation at the moment t. 

Here one can take as sq the critical size and the time to has to correspond 
to the time of the end of relaxation (actually it is the time of relaxation). 

To observe the boundary condition f(s — 0) = it is necessary to take 
the combination 

fo = /+ + /-, f+ = G(s,t\s ,t ), f- = G(s,t\ - s , t ) 

This gives the solution of this problem. 

Consider the behavior of the critical size p c . It is important to know p c 
because u = p/p c . One can propose the equilibrium critical size p ce which 
can be determined on the base of the size spectrum as 

r f ( P ,t)3p 2 (i-^)dp=o 

J— oo p 

This corresponds the stationary value of the critical size, i.e. dp c /dt = 0. 

It is clear that never p c equals p ce because this means the stationary 
value of p c and of the supersaturation (. But approximation p c p ce is 
rather good. Namely, at the beginning of the diffusion blurring this equality 
takes place. So, the critical radius p c (we shall mark it p c0 for initial time) is 
given by condition 

poo 

/ fo(p)v p (pco)dr = 

J 
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Here v p is the velocity of growth of the variable p. It is supposed that the 
size spectrum is relatively narrow. Another variant taking into account the 
different volumes of embryos is the following 



and see that ordinary y is small. Then one can analyze the evolution of the 
system through decomposition on y. 

When y is big it means that the size spectrum is essentially greater than 
p c . But this corresponds to the evolution via supercritical embryos where 
we have extremely simple law of growth dp/dt = CAi- So, the combination 
of the consideration of supercritical embryos and decompositions on y with 
restriction in the first several terms (actually the first non-zero term) can be 
very effective. 

Now we return to the diffusion blurring. The result of the diffusion blur- 
ring is very optimistic for the final conclusions. It sounds as following: The 
distribution function is the universal one and does not undergo the change 
of the form any more. So, one can see that the asymptotic solution is found. 
But the situation is not so simple. 

Really the function f after scaling in units of p ce will be the universal 
function without any parameters. 

On the base of distribution we can calculate the behavior of the critical 
radius. Note that it is impossible to find the critical size directly from the bal- 
ance equation 2a/3p c + / °° p 3 f(p, t)dp = const because here p 3 f(p, t)dp ~ 
const and the error radically increases. Instead of the direct balance equation 
one can use the differentiated variant 2(a/3pl)(dp c /dt) = 3 J °° p 2 f(p,t)dp 
which allows to find p c 

If the size spectrum will be the universal function then the critical size 
will be also the universal function. 

The special question is the correct boundary condition at small sizes. 
The velocity of the dissolution of the small embryos is a rather complex 
function of size and the regular dissolution exists. One can not neglect this 
regular dissolution. But fortunately the small embryos are dissolved very 
quickly with the growing velocity. So, one can suppose that they disappear 




Also one can propose to extract the deviation of p c from p, 



ce 



V = (Pc-Pce)/P, 



'ce 
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immediately at p — pf — (0.6 -r- 0.8)p c . So, the zero boundary condition has 
to be put at pf and all further considerations remain without reconsideration. 
Certainly, to keep the boundary condition we have to put the negative Green 
function symmetrical to s with respect to Sf. 

In the situation with Sf we have also the universal spectrum and the 
universal behavior of the critical size. 

But this universal asymptotics is only the intermediate asymptotics. 

Now we shall introduce the regular growth and destroy this universality. 

It is necessary to put some boundary of the type 



and for the sizes greater than p r one has to consider the regular motion with 
the asymptotic velocity. The choice of p r can be made also on the base of 
the LS analysis. 

It is trivial to refine the solution by investigation of the transition zone 
explicitly. 

The growth of the supercritical tail leads to the growth of the size p c which 
can be calculated in approximation of the following iteration procedure: On 
the base of initial p c0 we find the tail j ta ii of the size spectrum 



(the explicit dependence of v p on t is weak or it can be expressed via t 
iteratively) . 

The new distribution function f\ will be the superposition of the initial 
part /o and the tail. Then from the balance equation 



we find p ce i. This will be the new equilibrium critical size. It will be near the 
real critical size unless the tail begins to play the main role in the metastable 
phase consumption. 

The period of the diffusion blurring come to the end when the velocity 
of the growth for the critical embryo becomes to be comparable with the 
velocity of the growth for the tail. 



p r = (2 4- 3)p c 



ftail(P,t) = fo(pr,t')v p (p r )/v p (p) 



where 
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One has to analyze an attractive possibility to consider the process of 
diffusion in the region s < s r with the linear on s — s c rate of growth. Here s r 
has to be determined as the boundary between the linear and the asymptotic 
rates of the embryos growth. It seems from the first point of view that the 
linear rate of growth in the near critical region is preferable instead of the 
absence of the regular growth considered above. 

The Fokker-Planck equation under the linear rate of growth has the fol- 
lowing form 

dp = dyp | D & 2 p 
dt ^ dy dy 2 

Here the zero value of the coordinate as the critical size is taken, 7 and D 
are some constants. Here 7 is negative. The Green function for the positive 
7 at the infinite interval is well known and has the following form 

{X ' tlX ' ] ~ V 2tt£>(1 - e-27(*-f )) 6XP( 2D{1 - e-***-*')) ] 

Now it is necessary to take the combination of two Green functions and the 
answer is ready. The further analysis is absolutely the same. 

This approach seems to be more precise than the previous one but it has 
many disadvantages. The first disadvantage is the following: one can see 
that here 7 has to be negative and then at some time the half width of the 
gaussian goes to infinity. So, the solution becomes illegal. 

The second disadvantage is how to take into account the drift of the 
critical size. Now it appears in the rate of growth and then in the final 
formulas. The solution with a moving critical size is illegal also. 

But the idea to consider the law of growth as a combination of the linear 
dependence at p < p r and the asymptotic law at p>p r is very attractive. 
Really the term 1 — u~ x in the traditional law of growth can be treated as 
a correction term in the asymptotic decomposition. Here this asymptotic 
is taken over the positive powers of p r . The decomposition on the positive 
powers is not less justified in comparison with the decomposition on inverse 
powers. But the last approximate rate of growth allows an explicit integration 
and then the LS technique can be analyzed explicitly. 

One can try to construct the approximate Green function for the case of 
the presence of the regular growth in the following manner. We construct 
this function for initial perturbation at x = which is an equilibrium value 
for the regular growth v x (x = 0) = 0. We suppose that there are no other 



17 



points where v x = and v x (x) = —v x (—x). Then we seek the Green function 
in the ordinary form 



G(x, t\to) = A(t, t ) exp( 



x 2 



A(t,t ) 2 ' 



where the amplitude A can be reconstructed on the base of normalizing 
equation / Gdx = 1 and the width A is found by relation 



A = yjAD{t-t )+ t v x {A(t'))dt' 



or 



A = D(A(t'))dt' + J v x (A(t'))dt' 



'to •'U) 

for varying D. The last equation can be easily solved iteratively 



A i+1 = J A f DiA^dV + f v^A^dt' 

V J to J to 

A = A(t ) 

To ensure the correct boundary condition it is necessary to add the sym- 
metrical negative Green function. 

To refine the solution one can also use here the values of effective diffusion 
coefficient and effective law of the regular growth velocity from consideration 
made in [9]. 

4.3 The period of the dissolution of the head of the 
size spectrum 

The tail of the size spectrum grows and earlier or later the main role of the 
metastable phase consumption will belong to the tail. This opens the period 
of dissolution of the head of the spectrum. 

This period allows a rather trivial description since the high accuracy is 
not important here. Inevitably the head of the spectrum will be dissolved 
and this marks the end of this period. 

The most primitive description is the following. We split the substance 
between the tail G ta u and the head Gh ea d- The spectrum in the head is 
described by f . The spectrum in the tail f ta u is the direct translation of 
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the blurring part of the head which comes to the zone p > p r by the regular 
growth with the known supersaturation. The supersaturation corresponds in 
the first iteration loop to the stationary position of the critical size. Then we 
can calculate the dissolution of the head, the growth of the tail and replace 
the critical size on the base of the balance equation. We know a new value 
of supersaturation. This closes the iteration loop. 

Here it is impossible to use the model of the growth with a zero value 
of growth for p > pf and the zero value of p for all p < pj (otherwise 
it produces the jumps in the supersaturation value). One has to use the 
explicit law dp/dt = (2a/3p c )(l — (p c /p))- 

Another style of description is to use the methods from description of the 
dissolution of the tail of the size spectrum which is analyzed below. Since 
the method is the same we do not consider it here explicitly. Certainly, the 
exponential tail like exp(— const p) has to be changed to the head of the size 
spectrum f . 

This period is rather short and it is not very important for the further 
evolution. Details of this period description can be found in [8]. But one can 
see that in the theory presented here the new head at the tail is not formed. 
Here lies the main difference between this theory and the theory from [8j. 
The question whether the new head is formed is rather complex. This is the 
question of applicability of the gaussian tails of the Green function of diffusion 
equation. If we adopt the model with a finite upper limit of the size spectrum 
then we have a new head at the tail and have to use scenarios proposed in 
[8J . If we believe in long gaussian tails we come to the theory presented here. 
To solve this question we must go ahead of the level of description taken 
in the diffusion approximation. Otherwise there is no sufficient statistics 
to solve this question. We prefer to stop here at the statement that there 
is no sufficient statistics. It means that concrete details will determine the 
scenario. For example, we adopt that the act of the molecule consumption by 
the embryo requires a certain elementary time and we come to the finite upper 
limit of the spectrum. In the opposite situation we come to the gaussian tail. 
Certainly this question is out of the level of consideration adopted in the 
nucleation theory. 

4.4 The period of the gradual consumption of the tail 

Now we come to the period which in some sense resembles the LS results. 
The result of the previous period is the formation of the exponential tail 
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at p > p r . Now we consider the process of the tail dissolution. Here one can 
see the true competition between the embryos with different sizes. So, here 
it is convenient to go to the LS coordinates. 

To investigate this period we simplify the rate of the embryos growth. 
We assume the following rate of growth 

du _ ± 

— ~ (1 — u ) — 7« 

GST 

If there is a sufficient tail and the heat is already dissolved then at the 
moment of the end of the previous period 

f du 

max{ — } = v m < 
dr 

So, we assume that v m is negative. 

We split the whole region of u into three small regions. In the region of 
big u we suppose 

du 

— ~ (1 — 0) — 7« 
dr 

This law allows integration even under the variation of 7. 

The next region is the region of the intermediate u. Here we assume 



du 



= V m 



In the region of small u we neglect and get 

du _! 
Tt ~{1-u ) 

This law does not contain parameters and can be easily integrated. It means 
that the dissolution here is so fast that we can neglect the change of the 
critical size. 

The boundaries ui and u<i between these regions can be established from 
the continuity of the rate of growth. 

One can also put an effective boundary of the total dissolution instead of 
the zero size. 

It is possibly to refine the law of growth having introduced instead of v m 
some other effective value of the flat region. 
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When the size spectrum is known then the balance of the substance be- 
comes the transcendental equation on 7. After we found 7 we can get p c by 
integration. 

One can easily follow the dissolution of the spectrum on the base of the 
approximate rate of growth. 

We accumulate the approximate law of growth in the following formula 



du/dr ~ [du/dr] 



appr' 



This law of of growth allows to know u[r) on the base of some u [r ) for 
every arbitrary curve 7(V) 

u(r) = u (r ) + f [du/dr) appr = F appr [r\u , r ; 7(7-' )) 

For the law of growth we write 

du/dr = ip(u, 7(7")) 

Here we can take both approximate or precise law of growth. 
The substance balance equation 

dp c 2a 



roo 

dr 3p2 (7y = Pc( T ) J Q 3F a PP r( T \ u o, r ] d In Pc (T')/dT') 

^{F a ppr{r\u , r ; dhi p c {r') / dr'), din p c /dr)f (u , r )du 
dp c f°° 

+ 3 Pc( r )-^T J Q F appr( T \ U 0,To;d\np c (T')/dT')fo(u ,T )dUo 

is now the closed equation on p c (r). Here instead of a one can put the 
appropriate constant in accordance of normalization of the size spectrum. 

It is necessary to stress that all functional dependencies here are explicit 
ones and except p c (j) all other dependencies are known. The best way to 
solve this equation is to use the steepest descent method. The methods to 
solve this equation will be discussed below. 

This equation can be approximately simplified. Since the last term of the 
rhs is many times greater than the lhs one can approximately write 

POO 

Pc(r) J o 3Fa Ppr (T\u ,T ;d\np c (T')/dT') 
^(F appr (T\uo,r ;d\np c (T')/dr'),d\np c /dT)fo(u ,r )du = 

~ 3 ~dr Jo ^pp^^ ' T o'^ dln Pc( T ')/ dT ')fo( u o, r )du 
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or 



Jo^^appr^K^o; d\np c (T')/dT')(p(F appr (T\u ,T ]d\np c (T')/dT'),d\np c /dT)f (u ,T )du 
3 io°° F a PP r( T \ u o, r ;d\np c (T')/dr')f (u , T )du 

d In p c 



This equation can be solved by the same methods but it is more simple 
than the previous one. We outline again that except j(t) all other depen- 
dencies here are known. 

Another possible approximate variant of the balance equation is the fol- 
lowing 

PCX) 

Pc 3 = const' 1 F%pr{T\uo, T ;d\np c (T')/dT')f (u , T )du 

It seems to be the most simple variant of the balance equation. 

Now we shall discuss the asymptotic properties of the balance equation. 

Generally speaking the problem is solved since we know the good approx- 
imation 

d In p c /dt = 7 + some positive small value. 

We can solve this equation by decomposition in series or by some effective 
linearizations. 

But below we shall analyze the properties of solution in order to see that 
the size spectrum here resembles the LS theory for the size distribution. 

First of all we have to note that the tail has the exponential character. 
Really, the translation of the gaussian at some shift from the maximum leads 
to 

ftaii ~ exp(-constft) 

which can be easily approximated by the standard exponent of the argument 
linear on the size. Here the const is some fixed value proportional to (s r — s c ) 2 . 

Here we have to recall that in the original paper by Lifshic and Slezov 
there is a reference on the exponential tail of the size spectrum. It is quite 
natural to check the theory on example of the exponential tail. 

The exponential on r spectrum is exponential on u also in the asymptotic 
limit. 

In frame of the steepest descent methods the utilization of the exponential 
tails is quite justified. One can simply refer to the steepest descent method 
instead of the explicit consideration made above. But one has to stress 
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that we follow the explicit determination of the size spectrum instead of the 
formal methods. Explicit decompositions also give the exponential tail of the 
spectrum. 

If the characteristic width of the tail is many times greater than the 
critical size then v m is far from zero and there appears the rapid dissolution 
of the size spectrum. If the width of the spectrum is many times less than 
the critical size then v m is close to zero. It is evident that earlier or later 
the last situation will take place. One can give the qualitative picture of the 
process - The evolution at the big finite time is the slow monotonous increase 
of v m up to zero. 

The situation of the wide tail can be investigated rather elementary. The 
behavior of supersaturation is governed by the consumption of the substance 
by the wide tail. To see this consumption one can use the regular growth. 
This is described in [8] under the investigation of the oscillating regime. 
Evidently, the consumption of the substance leads to the growth of the critical 
embryo and the dissolution of the tail. This process will take place until the 
tail (or the rest of the tail) can be considered as the wide one. 

The rest of the tail earlier or later will become the narrow tail and then 
we can use the theory of the narrow tail. 

Now we consider the situation of the narrow tail. 

We consider the form of the size spectrum. We use the known formula 



where v u is the velocity in u-axis, t(u) is the time for the embryo to attain 
u. 

One can note that 



This asymptotics makes the size u inconvenient for analysis. It is more 
convenient to act in the p-scale where the asymptotic rate of growth is the 
constant one. In experiment under the instantaneous observation the variable 
u is proportional to p and there is practically no difference between them. 

So, it is preferable to consider at big u the p-scale. In the variable p 
the picture is rather simple - the exponential tail begins to be transformed 
according to the variation in the velocity of growth 




f(p) 



exp(-constp) 



(1 " Pc/P) 
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The amplitude of the spectrum and the width are determined from the 
behavior of p c (t)- Because of the asymptotic neighborhood of the behavior 
of the critical size to its behavior is the LS model these characteristics are 
close to the results of the LS technique. 

Now we turn to the justification of the similarity of the form of the size 
spectrum in the LS technique and the spectrum established in this theory. 
Until u ~ u m or p ~ u m p c there is no spectrum in the LS theory. In the 
current theory the situation is analogous - the tail is very short. 

Now we investigate the region u pa u m . Since u m is big the rate of growth 
dp/dt is close to the asymptotic value, i.e. to the constant velocity and the 
tail in the current model will be close to the exponential one. But what will 
be in the LS theory? We turn to the formula 

d(r - t(u)) = #(exp(r - r(u))) 

which can be rewritten with account of an initial form of the size spectrum 

as 

6{t — t(u)) = exp(const(r — t(u))) 

Having recall that 

. . f u du 

t(u) = — h const 

Jo v u 

which gives under the constant value of v u the evident relation near the 
maximum 

t{u) ~ u + const 

we see that the dependence of 6 on u (r is fixed but it is excluded) becomes 
the exponential one. 

Certainly here the derivation differs from the LS analysis and we ignore 
the change of v m in time which can be very essential. But qualitatively we 
come to the same results. 

The distributions in the region with small u are formed both in the LS 
theory and here by the dissolution of the exponential spectrum. They are, 
hence, similar. 

When one neglects 7« in comparison with (1 — u^ 1 ) it means that we 
neglect the change of the height of the original spectrum because the time 
of the dissolution of the given embryo from the size u ~ u m is small and the 
change of the critical size during this time is small. This simplification is 
quite possible. 
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So, both distributions (in the LS theory and the derived here) are similar 
in their dorm. The similarity is ensured by the narrowness of the tail of the 
size spectrum. Namely the situation of the narrow tail is the dominating 
one in the evolution scenario. So, the similarity is the occasional coincidence 
corresponding to the initial exponential form of the size spectrum tail. But 
namely this coincidence leads to the experimental confirmation of the LS 
technique. 

Later we return to the situation of the wide tail. Every wide tail as the 
narrow tail is also local in the size axis and, hence, there is the backlash 
in dependence of v u on time. The tail can be approximated in frames of 
the steepest descent method by an exponent. So, the style to construct the 
solution will remain the same. Hence, everything is reduced to the already 
analyzed situation. 

One can see the following stabilizing property - the wider is the tail, the 
wider is the backlash and the local character is approximately conserved. 
This property is very important - it is responsible for the observation of the 
LS-like spectrum already at the moderate time. 

Here the free molecular regime is adopted, this allows to write the balance 
equation in the integral form. The opposite case is the case of the profiles of 
metastability around the embryos. In this case one has to take into account 
the interactions between these profiles. The task seems to be extremely 
complex. Nevertheless the answer for the form of the size distribution is very 
simple. Certainly this answer is rather approximate. 

Really, one can propose the following model. Since the profiles are sharp 
functions of the space coordinate one can imagine only the pair interactions. 
Such a pair battle will end by defeat on one of partners. The winner will 
continue to be the embryo of a new phase, the looser disappears. At the 
asymptotics of evolution the remaining embryo had to win many battles. 
With probability pi it wins the first battle, with probability pi it wins the 
i-th battle. The total probability P to t to win all battles is the product YiiPi 
of all probabilities. Since pi are independent stochastic values we have for 
P tot the log-normal distribution. 

Certainly this approach can be spread to the group interactions (triple, 
etc.) which will give the same final result. 
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5 Development of the model 



Now we can note several important properties of the size spectrum. The first 
property concerns the influence of the boundary condition on the tail of the 
size spectrum. Consider 

/=-/- + /+ 
/_ ~exp{-{x + x ) 2 /4Dt) 
/+ ~cxp(-(x-x ) 2 /4Dt) 

Then 

/_ ~ exp(-x 2 /ADt) exp(-2xx /ADt)exp(-x 2 /ADt) 

f+ ~ exp(-x 2 /ADt) exp(2xx /ADt) exp(-x 2 /ADt) 

/+//_ ~ exp(4xx /4Dt) 

and one can take into account in the tail only the term /+. 

One can add that the regular growth can not destroy the tail - one can 
speak only about the shift of the tail and the sequential cut-off of the regions 
preceding the tail. 

The second property is the possibility to sweep out the boundaries be- 
tween stages in the sequential description of the evolution. Really, does the 
diffusion stop after the end of diffusional blurring? Certainly, it continues and 
the time of diffusional blurring depends on the amplitude of the spectrum, 
i.e. of the rest of the tail. 

One can note the following important feature - The tail blurring is so fast 
that it can not be overcome by the regular growth. So, the time of diffusion 
blurring is important and the diffusion process occurs during the whole time 
of evolution. 

Now we specify the recipe of calculations for concrete case. We shall ex- 
plicitly see what effect has the relatively small backlash in the law of growth. 

5.1 Explicit calculations 

We start from the law of growth 

^ = I(l-„-i) 
dt ti 
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for u = p/pc- Since 

dp dp c u du dp c 
~r — — ; — = Pc~, — V u—r- 
dt dt 1 dt dt 

du _ J__C _ _ u dp c 

dt p c t\ p c dt 

_ 2a 

du 2a . u dp c 

~dt ~ ^ _M ' ~p c ~dJ 

3f%h du _ . _ _ 3pchdp, 
We introduce r to have 



we see that 



Since 



we come to 



or 



u 

2a dt ' 2a dt 



2a , 

dt = dr 



3pl(t)t 



and then 



or 



du . 1 dp c 

T = ( 1 -^ 1 )-—T u 
dr p c dr 



So, 
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Now we find the argument u m which provides maximum for the rate 
of growth du/dr, i.e. the maximum of the curve (1 — u^ 1 ) — ^u. Having 
differentiated du/dr on u we have 

— [(1 — u^) — 7u] = u~ 2 — 7 
aw 

Then 

u m = T 1/2 

The height of the curve (1 — u^ 1 ) — ^u will be 

du. ■, 

— i _ 9?/ 
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It has to be zero or some small negative value —5. Namely 5 is the backlash. 
Then 

1 - 2U' 1 = -5 

and 

2 



7 = 



1 + 5 
(l + o) 2 



4 

Now we reconstruct the dependence of p c on t based on the known value 
of 7. We have 

(l + o") 2 d\np c 



or 



Then 



4 dr 

3tidp*_ _ (1 + 5) 2 
4a dt 4 



1 w^ + s ^ h 



'c 



or 



Then the supersaturation satisfies the asymptotic behavior 



Aati 



I 3(1 + 5)H 

5.2 Contradiction in asymptotics 

Now we can see the concrete picture for the approximate law of growth. 
The asymptotics 1 — 7« at big u crosses the axis du/dr = at 

4 

u — 7 



(l + o) 2 

The level of the backlash —5 it crosses at 

«r = (i + aft" 1 = Tj^y 
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The asymptotics I — u 1 for the rate of u growth at small u crosses the 
axis at 

u = 1 

and crosses the level of the backlash at 

^ = (l + 5)- 1 

So, the approximate rate of growth is constructed. One has to take into 
account that the backlash —5 can be also the function of time t or r. 
Now we analyze the asymptotic behavior of big u We have 

du 

— ~ 1 — ju ~ — ju 
dr 

Then 

lnw ~ — 7T u ~ exp(— jt) 

Now we can explicitly express r on t based on the known asymptotics p c . 
Really, 

2a , 

dt = dr 



or 



2a , 

dt = dr 



3(ft)*i 



Then 



and 



2 , 

-at = dr 
t 



2 In t ~ r 

Then the asymptotics for u in t— scale will be 

u ~ exp(— 7T) ~ t~ 27 ~ t -1 ^ 2 

Now we can see what will be the diffusional blurring exp(— const s 2 /t). 
We see that 

exp(— const s 2 /t) ~ exp(— const p 2 c u 2 / 't) ~ exp(— const -^—tt~ l /i) 

ot\ 

and it seems that the diffusion blurring is the main effect. It is no more 
than an error. The reason is that s does not grow here. Really, p = p c u ~ 
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£i/2£-i/2 _ cons i does n0 ^ g row . This occurs because we throw away the 
negligible constant in the law du/dr = 1 — ^u. But the effect of non-zero 
growth of p manifests in the constant 1 in this law. Certainly, it will be lost 
in comparison with the leading term. So, we see that the variables in LS 
theory are very dangerous. It is forbidden to choose new irregular variables 
and then fulfill the asymptotic analysis. 

As for the smallness of diffusional blurring one can easily see it directly. 
Since dp/dt = CAi an d we already know that £ ~ t^ 1 ^ 2 then the integration 
gives p ~ t 1 / 2 and s 2 ~ t 2 ^> t. This shows the smallness of diffusional 
blurring. 

One can take into account the modifications of the model 

• At the tail of the diffusion gaussian the regular shift is not very impor- 



• The initial diffusional blurring does not stop at the beginning of the 
regular dissolution of the tail but takes place all time long. 

In the LS approach it is used that all substance is in the region less 
than u m without justifications. Here we shall show the analogous fact (all 
substance is near u m ) explicitly. This fact has to correspond to the smallness 
of the diffusion blurring. We have to show this smallness. Really, the width 
s 2 ~ t of the diffusion blurring is many times less than the critical size 
p c ~ t 1 / 2 since s ~ p 2 . 

5.3 Modifications of the model 

The balance equation is the main instrument to determine the evolution of 
the system. It can be written in the following form 



Precisely speaking one has to add the the supersaturation as 2a/3p c and get 



but the first term goes to zero. In any case it is impossible to determine 
the critical size from the last equation having calculated the integral term in 
some approximation. 



tant. 
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We see that in analysis of the balance equation lies a dangerous pos- 
sibility to get wrong results. This possibility is extremely high in the LS 
analysis where the size spectrum has to cancel divergence of the critical size 
asymptotic behavior. 

To see the behavior of p c we have to establish the form of at least in 
the asymptotics. We write 

4>{U,T) ->■ </>aa(u,t) 

where as is defined at big u from the gaussian (one can show that the back 
side gaussian is not important here). Namely, we have the following chain of 
equalities 

(p as (u,T)du = f(p,r)dp 
f(p,r)dp = ^(s,r)ds 

* (s ' T) ~ eXp( -4^ } 

where s is p 2 , D s is the diffusion coefficient over s (known value). 

Roughly speaking the problem is solved. But we can not combine the 
values of u, p c and r until we integrate the law du/dr = (1 — u~ v ) — 
of growth. Fortunately we can not do this with 7 varying in time. So, it 
is necessary to introduce approximations for this law of growth. This has 
been done above. Actually we are interested now in the size spectrum for 
the values u > u r . Then we have 

r(t) du. , rr<t). 
u = u+ ( — )dT = u+ / (1 — 'jujdr 

Jr(t ) dT Jr(t ) 

where u is the size of uaXtQ. Here to is the time when the diffusion transforms 
into the regular motion and u is the corresponding size. The value of u can 
be found on the base of u, r and the initial value, then we take the initial 
size spectrum at u and get the size spectrum for u at r. 
Very approximately we can substitute the law 

du 

— = 1 — 7« 



by 

du 
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This gives 

r r(t) 



rT(t) 

■u = wexp(— / 7(r)dr) 

Jr(to) 



or 

/•r(t) 

u = u exp( / j(r)dr) 

Jr(t Q ) 

But as we have seen earlier this leads to an error. Fortunately we can inte- 
grate already du/dr = 1 — ^(t)u explicitly without simplification (formulas 
will be long). 

On the base of u one can find 

P = Up c (t ) 

and then s equal to p 2 . 

For the distribution in p-scale it is very easy to write the gaussian 

n i P 4 \ d P / P 4 \ 

G ~ exp — - — ~ exp — — — — ) 

PV D s t(r) ) ds PV D a t(r) J 

Here we ignore the jacobians arrived from transition from s to p scale because 
exp(— constx A ) at the tail is a very sharp function. We shall ignore them 
below also. 

We are interested now in behavior at u ~ u r . Here the true approximation 
will be 

du 

Tt =C4 ~ lu 

Ci = 1 — u^ 1 

Let us take 7 ~ 7^ = 7(^s) (^b is the time when u r is attained). 
Then 

u - ci/7b 

So, 

« - Ci/tb = (« - Ci/7b) exp(-7 B (r(t) - r(t ))) 

and -u is u at t - As the result we know u as a function of w, i.e. u = F(u). 

In any case we can integrate the equation with 7 = 7(7-) since we have 
the first order linear differential equation. This will give the real true result 
with rather long formulas. So, we use the previous formulas keeping in mind 
the necessity to apply the formulas with varying 7. 
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Now we can see the asymptotics for the distribution function 

u{uYp A c ci - -f B u(u) 

<P(u,t) ~ exp 777^ — 

D s t ci - 7 B m 

This result is for u > u r . 

One of important properties of this solution is the weak dependence of 
the form (after scaling) of the size spectrum on r. The quantity of substance 
in the region [u r , 00) can be found as 



G> = pi u 3 (f)(u,T)du 

J U r 



Now we shall analyze the region [ui,u r ]. The values at this region are 
marked by the subscript =. Here 

du _ 
dr 

Then 

du r (f ) 

(f)(u,r) = 0=(m,t) = <j>{u r ,T)-£fa- 

where f is defined as r 

u = u- J 5(r')dT' 

which is a rather complex closure. The previous equation can be rewritten 

as 

<f)( U ,T) = 0=(m,t) = 0(w r ,r)— - 

6{T) 

More rigorous is to make a shift of u r to exclude a slow evolution right of 
u r . But this does not change the qualitative behavior. 
The quantity of substance G= in this region is given by 



/ttj ./It; 



Here 

■u r = u + 



f T 5{T")dr" 
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Now we investigate the region [0,1^]. The values at this region will be 
marked by the subscript <. Here the size spectrum can be given by 



^(«,r) = ^(u,y) T ^ I 
1 — u 1 

The quantity of substance is given by 



du 



Here 

r u i du 



-r'=r- 

Ju 1 



u- 1 



The unknown function is S(t). 

All quantities of substance G < , G=, G > have to be substituted into the 
balance equation which gives the equation on S(t) with the known coefficients. 



5.4 Steepest descent procedure 

Now we shall discuss the ways to solve this equation. At first we have to get 
a true algebraic equation. We differentiate the balance equation on time and 
get 

d(G < + G = + G > ) _ Q 
dr 

Then we use concrete approximations to get dG > /dr, dG < /dr, dG = /dr. 

For dG > /dr we see that the subintegral function is the product of the 
three functions: 

1. the moderate function 3u 2 , 

2. the rapidly growing function for the absolute value of the rate of growth 
du/dr. This function becomes very small at u r (practically it is zero). 

3. the rapidly decreasing function for the initial size spectrum. This func- 
tion decreases in the main term even faster than exp(— x 4 ) and even 
being multiplied by the accelerating rate of growth and by u 2 the prod- 
uct goes to zero at big u. 
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So, one can effectively use the steepest descent method. 

For dG = /dr we have a simple expression which is approximately propor- 
tional to the size spectrum multiplied on 3u 2 5 and due to the rapid decrease 
of the size spectrum we see that the subintegral function is the rapidly de- 
creasing function of u. This is the ideal situation for the application of the 
steepest descent method with the maximum at the boundary point (here it 
is ui). 

For dG < /dr we have the integral with the subintegral function at the 
interval ]0,W/], which is the product of three rapidly varying functions: 

1. The function 3m 2 which goes to zero at u — 0, 

2. The function du/dr ~ 1 — u^ 1 which goes to infinity at u = while 

u 2 (du/dr) goes to zero at u — 0, 

3. the rapidly decreasing size spectrum. 

This provides good conditions for the application of the mentioned steep- 
est descent method with the point of decomposition inside the interval. 

The problem which can appear here is a too sharp form of the size spec- 
trum which can cause the maximum of the subintegral function at p = up c 
which is too small. Fortunately when such sharp size spectrum will come 
to the region [0, u{\ the decrease of the substance in embryos will cause the 
decrease of p c . It means that until the diffusion blurring the spectrum will 
not be dissolved. This situation is typical for the formation of the new head 
which has been discussed in [8]. 

The mentioned approximations lead to the algebraic equation on 5. This 
equation allows the further simplification if we decompose 5(t) in Taylor's 
series on inverse time £ = r _1 

7(0 = 7o + E^e 

with coefficients aj. The choice of £ as a variable is recommended by the 
structure of the correction term in the LS procedure. Also one can use 
decomposition near some value £ 

7(0 = 7(£o) + £«;(£-&) j 
j 

The last modifications make the analysis of the evolution a technical task. 
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6 Situation of wide tails 



Now we shall analyze the situation of the wide tails of the size spectrum. It 
means that the half-width of the size spectrum in the region p > p r is many 
times greater than p c . It is clear that the approximate law of growth du/dr = 
1 — 7« corresponds to the law of growth dp/dt = (/ti and the distribution 
function moves as a whole along p-axis. The critical size p moves faster (but 
with the same time dependence t 1 / 2 ) and eats the spectrum sequentially. 
For all distribution tails (here we take the tail multiplied by the number of 
molecules in the droplet p 3 ) which decrease like 1/p 1 with % greater than 1 the 
relative half-width will decrease. But the situation with % < 1 do not ensure 
the finite substance, it is forbidden situation. So, earlier or later the tail will 
be narrow. One can speak, thus, only about the intermediate asymptotics. 

Since here the width remains greater than the critical radius one can speak 
about another asymptotic behavior, at least the intermediate asymptotic 
corresponding to the wide spectrum can be observed. 

Kinetics of the process here will differ from the LS case. All the time the 
main consumers of the metastable phase substance will be the big embryos 
and the backlash here is very wide. The asymptotic 

t(u) = 7 In w 

ensures the infinite time of the embryos dissolution for the infinitely big 
embryos. Then the main supposition in the LS theory fails. This argument 
states that to keep the balance of the substance it is impossible that some 
relative size will only grow and it is impossible that all relative sizes will 
decrease. But now it is possible to see here the situation where all sizes 
decrease rather intensively but the main consumers of the vapor are the big 
embryos with u many times greater than u m . In the zero approximation the 
evolution is very simple - the size spectrum is cut-off by the critical size (for 
the sizes less than the critical one the dissolution is so rapid that we can 
speak about the instantaneous dissolution and neglect this region). This cut 
off corresponds to the conservation of the substance in the system. 

Define that here the size spectrum fi on g(p) has the characteristic width 
A(fi ong ) determined in the integral way as 

/ florMdp = fZ^iflong) 
J 
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where fj™f£ is the amplitude of the size spectrum, or in the differential way 

as 

fZsiPm + A(f long )) = /££/exp(l) 

where p m is the argument for the amplitude value of spectrum. This width 
satisfies the strong inequality 

A(flong) > Pc 

The substance balance here is written as 

$ = C + q+ 

where $ is the supersaturation without the embryos formation, q + is the 
substance in the tail which can be calculated as 

/*oo /*oo 
q+= flong(p)p 3 dp « / flong(p)p 3 dp 

Then in the last integral we can take for fi ong the size spectrum fully deter- 
mined by the supercritical law of growth from the "initial" spectrum f (at 
the time t init ): 

flong(p,t) = (( (Unit) /((t))f (p(t, t in it)) 

where p(t, t init ) is determined from 

P= f C(t')/tidt' + p(t,t init ) 

It is more convenient to study the initial spectrum and to see how much 
this spectrum is cut off. So, we introduce the initial p-size variable s and 
write a balance equation 



$ = C(t) + / y?( S ) 3 /o(sK< 

J Sr. 



Here s c is the initial size of the variable s which attains at t the size (2-=-3)p c , 
<p(s) is the size which will be attained at t by the embryo with initial size s 



(p(s) = s+ I 
Jt, 



^ tinit ^1 

The balance equation can be rewritten as 



$ = C(*) + / ft c(t0 , , , 2H ^fj{s)ds 
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or since (2 + 3) has to be many times smaller than the width of the 
spectrum then 



t init '1 



Having noticed that s has to be at least positive we get the following modi- 
fication of the balance equation 



$ = at) + [ V (s) 3 f ( S )ds 

Js>0 



This equation can be easily solved. Having introduced the explicit equation 
for ip we come to 

* = <(t)+[ (s+ f ^p-dt') 3 f (s)ds 

Js>0 Jtinit II 

One can see that the integral term is the polynomial on 



Then 



with known constants 



It can be rewritten as 



Pm= ——dt 

Jtinit Cl 

* = «<) +t^f^A' 

a-i = / s l f (s)ds 
Jo 



The last equation is the differential Abel equation - the ordinary first oder 
differential equation. Since there is no explicit dependence on the argument 
this equation can easily integrated. This gives the solution of the problem. 

Certainly the last solution is not accurate because there all embryos 
remain in the integral term - they remain supercritical ones. This leads to 
the qualitatively wrong behavior. It is clear that the error will be essential 
namely when the tail stops to be really the wide tail. But nevertheless it is 
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possible to take into account the dissolution of the embryos. The balance 
equation has to be written as 

S = C(f)+ / (8+ f ^P-dt'f f ( S )ds 

and conserves the polynomial structure on p m . 



C(t)+ [ {S + Pm ) 3 f0{s)ds 

Js>(2+3)2a/3C 



or 



/s>(2H-3)2a/3C 

3! 



; j!(3 - Jj! A>(2-3)2a/3C 



Now the coefficients 



j!(3 - Jj! A>(2-3)2a/3C 

are known (since the initial size spectrum is known) functions of (. But it is 
possible to solve this equation on p m as the third power algebraic equation 

Pm = F(C) 

or 

p m = F{ti- 



dt ' 

with a known function F With the help of the inverse function F _1 we can 
write 

/.'/•• V.) 



dt 

The last equation can easily solved. 

Now we take into account the growth of the embryos. The balance equa- 
tion has to be written as 

* = C(*)+/ (*+/' ^f-dt') 3 f (s)ds 

Js+p m >(2^3)2a/3( Jt lnlt t\ 

and does not conserve the polynomial structure on p m . It can be written 
analogously 

* = C(fl+EA£ .,J* ! - M f , , , s3 - j f°( s ) ds 

j J -{3- J)\ Js+ Pm >(2-r3)2a/3C, 
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or in the standard form but with the coefficients 

3! 

j!(3 — j)\ Js+p m >(2~3)2a/3( 



J — ? ! Js+o m >(2+3)2a/3C 



depended (since the initial size spectrum is known) on p m and ( = t\dp m jdt. 
So, we have the first order differential equation without explicit dependence 
on the argument t 

^(p m ,dp m /dt) =0 

This equation can be integrated when we can express dp m /dt via p m . This is 
an algebraic problem which can be solved at least locally in a good approxi- 
mation. 



7 Approximations in the explicit construc- 
tion the size spectrum 

The first task in construction of the size spectrum is the construction of 
initial distribution which will be later gradually dissolved during the over- 
condensation. This task is very complex. Even the elementary approximate 
blocks for solution do not allow the solution. We do not know the solution 
of the diffusion equation for the diffusion blurring of the spectrum even with 
the stationary value of the supersaturation. As the result of such difficulties 
only very approximate methods can be formulated. 

To see the initial distribution it is more easy to use the p-scale because 
in this scale the asymptotic law of growth is rather simple and the velocity 
of growth contrary to the w-scale does not go to infinity. 

The first approximate model, which allows solution is "the model of se- 
quential evolution". We consider the period of the initial diffusion blurring. 
Here we have the stationary supersaturation (and the critical radius). 

Consider the value of distribution at some pf, time being fixed. Let the 
initial distribution be 5-function at p = p c . The route from p c to pf will be 
split between the diffusion blurring and the regular growth. In this model 
the blurring occurs up to p equal to some parameter pj,. Later there will be 
the pure regular growth. Parameter pb is reasonable to put equal to 2p c . 

The time of the regular growth up to pf will be (p/ — Pb)ti/( or (pj — 
Pb)ti3p c /2a. Then the time to quit the diffusion blurring will be 

t q = t- (p f - p b )ti3p c /2a 
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The supercritical regular growth is the simple translation of the size spec- 
trum, then we have to calculate the spectrum at pb and t q after the pure 
diffusion blurring (here for simplicity we do not consider the back side input 

/-) 

a, ,\ 2 Pt> / Pb \ 2 P?> / Pb \ 



yflDj KV AD s t q ' y/Wj KV 4D s {t- {p f - p b )t 1 3p c /2a) J 

The function d can be considered as the initial size spectrum. The prob- 
lem is, thus, solved. 

But the function d has a certain disadvantage - the size spectrum is finite. 
Really, for pf greater than pu m 

Plim = Pb + 2at/ti3p c 

the spectrum is zero. 

Actually the values of p near pu m begin to dissolve (i.e. pu m is near p r ) 
when the amplitude of the rest of the spectrum is extremely small. So, the 
relative quantity of droplets in negligibly small. It will be big only in the 
systems of cosmological sizes. 

Nevertheless one can refine this solution. Fortunately, the diffusion pro- 
cess can be easily estimated even with the varying p c because here ( is 
very small and the average number of collisions in the time unit will be 
(2/ti)(3p 2 (t)). Here we can take p on the base of the regular growth and get 
the total number of collisions 

n to t= l\p 2 {t')dt' /h 
Jo 

The characteristic width A tot can be estimated as A tot = y/2n to t with a 
sufficient accuracy. 

Then the resulting distribution will be proportional to 



= r eM^A^MPf^Pf 

J— CO ^lltnt 



Hot 

or 

^ = 2p b roo ^ (p~P/) 2 ex , pt )d 

y/mtJ-oo XPl 2n tot j6XPl 4 J D s (t-(p / -p 6 )t 1 3p c /2a) j Pf 
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It is quite satisfactory here to consider only the right hand wing of the dis- 
tribution, i.e. to put p > pj. Since the subintegral function is the product of 
exponents it is reasonable to use the steepest descent method. 

The required result is attained by a simple combination of the regular 
growth and the pure diffusion. The cross effects when, for example, the 
stochastic increase of the embryos size leads to increase of the regular rate 
of growth are missed here. They can be included into consideration by con- 
sideration of the effective half-width of the diffusion blurring and the shift in 
regular growth proposed in [9]. Certainly, they have to be slightly reconsid- 
ered since now we start not from the very beginning of the size axis. 

Now we shall show the primitive approximate way to construct the explicit 
form of the size spectrum for the narrow tail. 

The initial distribution is supposed to be known. 

We approximate the rate of growth for u in rescaled time r by the fol- 
lowing approximation 



du du 

— K,v m = max{ — \ = -d(r) 

dr dr 



for u G [ui, u r 



for u > u r 



du 
dr 



"fU 



for u < U[. 

Parameters u\ and u r have to be chosen to ensure the continuity of the 
whole approximation. 

Consider u > u r . The law of growth for p corresponding to du/dr ~ 
(1— 7«) is dp /dt = Q/ti (the r.h.s. is precisely taken into account by transition 
from t to r). So, the distribution in p-scale is moving as a whole without 
changing of the form. We know this form - it is exponential form 

f(p) ~ exp(-ap) 

with some parameter a. 

Then the distribution <p over u is connected with / by the following 
relation 

¥>(«) = fiP)t 

dt 
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Having inserted the explicit relations for the derivatives we come to 

1 -u- 1 



<p(u) = f(p)- 



(1 — u l ) — 7« 
or in the supercritical limit 

Here the distribution function is even sharper than the distribution over 

P- 

Now at first we shall assume that the backlash is changing in time slowly 
in comparison with the time of dissolution of an embryo from u T to U[. This 
case will be at least the base approximation for iteration procedures to refine 
the solution. 

The known solution in the region u > u r leads to the known rate of 
appearance \I/&(t) at u — u r . This value is given by 

^6(r) = ip{u)\ u=Ur ( — \ u=Ur ) = (p{u)\ u=Ur v m 

CiT 

Here tp is the distribution over u. Now it can be established. One can easily 
show that ^b( T ) at -u = -u r is sharp decreasing function of r. 

The last function is the source at the left side of the central interval 
[ui,u r ]. Now it is possible to solution at this interval. This solution is very 
simple and it is given by 

<p(u,r) = * 6 (f)/^| u=Ur 

or 

The time f satisfies relation 



u-u r = J_ v m (t')dt' m v m (r - t) 



Hence, the solution in the central region is constructed. 

Generalization for the varying 5 is rather simple. The rate of appearance 
^(f) at u = u r is given by 

%if) = vi u \ : T')\u=uV m {.r) 
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Here ip is the distribution over u. Now it can be established. 

The last function is the source at the left side of the central interval 
[ui,u r ]. Now it is possible to solution at this interval. This solution is very 
simple and it is given by 



Now only the last region - the region of small u < ui, has to be investi- 
gated. Here the solution is also rather simple - it is the simple translation 
of the source from u r under the law of growth independent on 7. Here all 
constructions are analogous to the previous case but the transition over the 
central region has to be taken into account. 

The distribution <p> is given by 



<p(u,r) = *e,(f)/u m (r) 



and 



<p(u,r) = ^(w|f)| u=Ur ( f )U m (f)/u m (r) 



The time r satisfies relation 




T 



<p(u,t) ~ * a (f)/(du/dr) ~ *„(f)/(l -u- 1 ) 



where 



^a(-r) = <p(ui(f),f)v m (f) 
Here t is determined by the following way 




and 




For S = const one can simplify the last relation 




or 
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The integral can be easily taken analytically and we get the explicit expres- 
sion for t. 

If v m does not essentially depend on time one can simply solve all these 
equations. The form of the size spectrum is determined, only parameter v m 
(7, u r , ui depend on v m ) is unknown. The balance equation will be algebraic 
equation on v m and can be easily solved since we know at least the zero 
approximation (for example, one can take LS asymptotic). 

It is clear that the backlash is not a constant value, it changes in time. 
The effective way to investigate the situation of the varying backlash is to 
consider this variation small, then to decompose v m (r) in Taylor's series, to 
take few first derivatives and to fulfill the same program as was done above 
in the case of the constant value v m . 

8 Finite number of embryos 

Under the finite spectrum of sizes the LS asymptotic will be also destroyed. 
After the size of the greatest embryo attains the critical value and 7 goes 
to zero the diffusion blurring begins and we return to the section about the 
diffusion blurring and the evolution will be described by the same formulas as 
mentioned above. Here it is important that the diffusion term plays the main 
role in evolution. Certainly, diffusion begins to play essential role earlier than 
the biggest embryo attains the critical size. 

We see that the further scenario is built on the doubtful alternative 
whether the size spectrum finite or not. The diffusion process through the 
formula for the Green function gives the infinite size spectrum. But is every 
concrete system the spectrum is the finite one. The answer on this question 
also determines the asymptotic. The type of asymptotic is determined by 
the time of observation and the sizes of the system under the observation. 
Fortunately this question is artificial because other effects (the change of the 
regime of growth, the thermal effects, etc.) lead to the end of applicability 
of the chosen physical model. 

In every system the finite spectrum is the direct consequence of the finite 
number of embryos. So, we have to develop methods to describe the evolution 
with the finite number of embryos. 

Suppose we have few embryos in the system. Then the evolution is deter- 
mined by the laws of their regular growth (diffusion has also to be taken into 
account but in the manner of some stochastic adsorption and ejection of the 
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molecules by the embryo). The balance equation links these laws of growth 
in the closed system of equations. The number of these equations equals the 
number of embryos. When the number of embryos is less or equal to few 
hundreds it is preferable to solve these equations by computers explicitly. 
We rewrite the law of growth for a chosen embryo in the following form 

dpi = (_ ^ _ Pc^ 
dt ti pi 

_ 2a 
Pc= 31n(C+l) 

(index i marks embryos), where a is the rescaled surface tension (it is a 
constant) and ( is the supersaturation. Namely through the supersaturation 
one can link the laws of growth having written the balance equation 

i 

Here <3> is the initial value of supersaturation. 

One can easily solve the system of these equations at least approximately. 
The first method is very simple. At given initial sizes of embryos we find the 
value of the supersaturation. Then we reconsider the sizes having moved 
them according to the law of growth at the given supersaturation. The time 
interval has to be small. Then we recalculate the supersaturation having 
made one step of evolution. This method is the step by step procedure. 

The stochastic adsorption and ejection of the molecules by the embryo 
can be taken into account by the following simple procedure. We keep in the 
memory of computer all coordinates of embryos pi and for every embryo at 
the given supersaturation we determine the rate of adsorption as 

R+ = 3p 2 (( + l)At/t l 

Here At is an elementary time step. The rate of ejection will be 

R_ = R+- (dp/dt)At 

dp = C_, _ 2a_ 
dt ti 3Cp 
The rate of staying still will be 

Rq = 1 — R— — 
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We must choose At enough small to have Rq > even for the greatest 
embryo. 

Having put the point at the interval [0, 1] stochastically we determine 
what action we have to do. If this point belong to interval [0, we eject 
the molecule, i.e. we make a transition pi — > pi — 1. If this point belong to 
interval [R-, R- + R+] we accumulate the molecule, i.e. we make a transition 
Pi —> Pi+1- If this point belong to interval [R_+R + , 1] we keep the coordinate 
still, i.e. we make a transition p, — > p^. 

We have to repeat this action for every embryo. Then we recalculate the 
supersaturation 

i 

and fulfill the step in time. 

We see that these procedures can not give the analytical properties of 
evolution with the finite number of embryos. Hence, the problem of the 
system of several embryos exists and the effective solution is absent. 

Now we turn to the simplest case - the case of small quantity of embryos. 
The asymptotic of the process is the evident there will be the greatest embryo 
which will be the critical one. This embryo is in the effective potential well. 
Description here is analogous to the case of the several identical embryos. 

Really, the law of growth 

dp $ - p 3 2a 

~dt = ty ( ~~ 31n($-p 3 + l)p' 

corresponds to the regular motion in the potential 

Jo h V 31n($ -p 3 + l)p ; 1 

which is very deep and only in the region of small p it has a barrier and 
begins to decrease when p decreases. 

If we ignore the fluctuational formation of the new embryos then the 
evolution leads to the fluctuational disappearance of the last embryo. It 
requires absolutely giant times and at these times new supercritical embryos 
have to appear. So, the ignorance of the fluctuational formation of new 
embryos (with extremely slow rate) is illegal. 

The flow of disappearance of the last embryo is many times less than 
the flow of appearance (formation) of the new embryo even in the system 
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with the practically consumed metastable phase by the giant embryo. So, 
one can see the formation of the second embryo. Later one can observe 
the competition between these embryos. Practically inevitably the second 
(new) embryo will be dissolved but with a small probability n which can 
be estimated (very roughly because we keep the old potential U appropriate 
only for one embryo) as 

7T ~ exp(C/(p = p min /2 1/3 ) -U(p min )) 

(here p min is the argument of the minimum of potential) the first embryo be- 
comes the second one and it will be dissolved. Certainly, the above presented 
picture is only the rough estimate. 

9 Conclusions 

In the theoretical constructions presented above we came to the following 
results 

• In Sections 2,3 we pointed out the weak features of the LS approach. 
In section 2 the features concerning the behavior of the supersaturation 
were outlined. In section 3 the weak features in construction of the size 
spectrum were presented. 

• In section 4 the sequential analysis of the evolution was given and it is 
shown why the form of the size spectrum resembles the form given in 
the LS approach 

• In section 5 the details of the most difficult and the most important 
period of the tail dissolution are given 

• In section 6 the approximate way to construct the initial size spec- 
trum for the period of the tail dissolution is given. This point is very 
important for qualitative results. 

• In section 7 the situation of the wide tail which allows essential simpli- 
fication is presented 

• Section 8 is devoted to the specifics of the case when only few embryos 
remain in the system. 
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Sections 7 and 8 are supplementary ones, the complete theory is given in 
sections 4-6. This theory gives an answer on two important questions - why 
LS spectrum of sizes can be really seen in nature and what is the difference 
between the real situation and the LS approach. 
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